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Abstract 

New physics contributions to B decays often arise through chromo-magnetic pen- 
guin operators. To look for new physics effects in B decays, it is useful to be able to 
estimate the hadronic matrix elements for the chromo-magnetic operator. We com- 
pute this contribution to B — > <f>K decays using PQCD methods. It is shown that, 
if the Wilson coefficient of the new physics is same order of magnitude as that of the 
Standard Model, this operator gives a non- negligible contribution compared to that of 
the Standard Model (about 30%). We also investigate the value of q 2 , which is the 
momentum transferred by the gluon in the chromo-magnetic penguin operator. We 
find that the expectation value (q 2 ) is approximately M|/4, in agreement with a naive 
guess. This result, however, is very sensitive to the scale dependence of the Wilson 
coefficient. We also show that the matrix element for the chromo-magnetic penguin 
operator is independent of the choice of energy scale to a very good approximation. 
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1. Introduction 



B — ► 4>Ks decay plays an important role in the search for new physics beyond the 
Standard Model (SM), as this decay is a pure penguin process. In the SM, the time dependent 
CP asymmetry for B — > (f>Ks decay is the same as that of B — > J/ ipKs decay, which is a tree 
dominant process. Any difference between them would be a signal for physics beyond the 
SM. Recently, the BaBar^ and Belle 2 ' collaborations reported CP asymmetry in B — > 4>Ks 
decay. If the present trend continues with additional data, the results suggest a deviation 
from the SM prediction. 

In the effective Hamiltonian approach^ the Hamiltonian is expressed as the convolution 
of local operators and the Wilson coefficients, which can include contribution from new 
physics. The effects of new physics change the Wilson coefficients. In many cases, the Wilson 
coefficient for the chromo-magnetic operator is most sensitive to new physicsP^^^^^^ 
The factorization approximation (FA) is often used in computing the hadron matrix elements. 
However, it is difficult to calculate this contribution in the FA, because the magnitude of q 2 
of the virtual gluon in the chromo-magnetic penguin operator is unknown. In the FA, q 2 is 
often assumed to be or Mfj/2. Then, the result is sensitive to the assumed value of 

In this paper, we calculate the chromo-magnetic penguin using the perturbative QCD 
(PQCD) approach within the SM. We can predict not only the factorizable contributions but 
also non-factorizable contributions and annihilation contributions, which cannot be calcu- 
lated in the FA. Furthermore, we can calculate strong phases, so that it is possible to predict 
CP asymmetries. PQCD has been applied to some exclusive decay modes, i.e. Kir^St 
KK, 14> Kr/, 151 D s KE® p(u)ttP (pR^BE® K*t^ and <j)K*f^ and the results are con- 
sistent with experimental data. They were calculated to leading order in a s in the PQCD 
approach. Therefore, the chromo-magnetic penguin was not included in the computation of 
the branching ratio for the B — ■> <f>K decay. Obviously, there are many other higher-order 
diagrams that must be considered simultaneously if we are to add the magnetic penguin term 
to the SM computation. To estimate the new physics contribution, the magnetic penguin 
contribution can be singled out, as it is most sensitive to new physics. In this study, we also 
analyze the FA computation of the color magnetic moment using PQCD as a guide. 

The outline of this paper is as follows. First, we present the effective Hamiltonian and 
the chromo-magnetic penguin operator. Then, we derive magnetic penguin contributions for 
hadronic two-body decays. Next, we calculate the chromo-magnetic penguin contributions 
in B — > <f)K decays using the PQCD approach and we give the numerical result. We also 
derive the distribution of q 2 . Finally, we summarize this study. 
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§2. Chromo-magnetic penguin operator 



In the effective Hamiltonian approach, the Hamiltonian is expressed as the convolution of 
local operators and the Wilson coefficients. The effective Hamiltonian for AS = 1 transitions 
is given by 



H, 



eff 



V2 



,q'=u,c 
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-V*V tb J2 CiMOiO*) + C 7l {n)0 7l {n) + C sG {n)0, G {n) 



,i=3 



h.c. 



(2-1) 



where V*, s and V q >b are the Cabibbo-Kobayashi-Maskawa matrix elementsp^^' Ox_io are 
local four-fermi operators, 7l is the photo-magnetic penguin operator, and 0§g is the 
chromo-magnetic penguin operator. The local operators are given by 





= {si<^) V -A(q'jbi)v-A , 


= 


= (Si<£) V -A(<fjbj)v-A , 


o 3 = 


(Sibj)y-A /X^jq^v-A , 
q 


4 = 


{sjbj)y-A yXq~jqi)v-A , 
q 


o 5 = 


(sibi) V -A^2(qjqj)v+A , 


o 6 = 


(sjb^y-A y^qjq^v+A , 
q 


o 7 = 


-{sik) V -A £ eqiqjq^v+A 
q 


, 8 = 


2 ^ 

-(sib^v^A^eqiqjq^v+A , 
q 


o 9 = 


q 


, o 10 = 


2 A 

-{sib^v-A^e^q-jq^v-A , 
q 


o 7l = 


g 


OsG~- 


= -|^m 6 ^(l+75)W^ 



(2-2) 



where i and j are color indices, and q is taken to be u, d, s and c. We use the leading 
logarithmic results of the Wilson coefficients in our calculations.^ We consider the chromo- 
magnetic penguin contribution for hadronic two-body decays. The non-local operator 0' 8G 
generated by 0$g is given by 



0' 8G = -^rn bSi a^T^{l + 75)^ 



—i 



+4ira s q»f abc [g u \q + k x )° + g x °{-k x + k 2 f + g™(-k 2 - q) x ] 



—1 
q 2 



x-rji (qvlxTi'j'qj') T2 i (qi'H<T T i!>j»qj") ( 2i 



(2-3) 
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where q is the momentum transferred by the gluon, which goes out from the a^ v vertex. The 
second line in Eq. ()2-3j) is induced by the self- interaction of gluons in G* , and the third and 
fourth lines contain a 3-point vertex of gluons. The quantities k\ and hi are the momenta 
that go out from the 3-point vertex of gluons. These terms are needed to maintain gauge 
invariance. The first term is first order in a s , and the others are of order a 2 s . In the PQCD 
approach, all of the terms contribute to the same order in a s , as we see later in the next 
section. 

§3. Magnetic penguin amplitudes in the PQCD approach 

In this section, we calculate the chromo-magnetic penguin amplitudes in B — > <fiK decays 
using the PQCD approach. The chromo-magnetic penguin operator for hadronic two-body 
decays is written as in Eq. ()2-3|) . and there must be at least one hard gluon emitted by 
the spectator quark in PQCD. Therefore, the diagrams with the chromo-magnetic penguin 
operator are as shown in Fig. ^ The diagrams (a) - (h) come from the first line in Eq. ()2-3|) . 




Fig. 1. Magnetic penguin diagrams in B — > cf)K decays. 

and the diagrams (i) and (j) come from the other terms. All of these diagrams contribute 
to the same order in a s . The chromo-magnetic penguin contributions are of next-to- leading 
order in a s within the PQCD formalism. However, if we include the Wilson coefficient, they 
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are of the same order in a s as the penguin contributions. If we consider the full diagrams 
of the penguin and the chromo-magnetic penguin, we see that they are of same order in 
a s . Therefore, we cannot say that the magnetic penguin contributions are smaller than 
the penguin contributions. In order to obtain a meaningful result for the branching ratios, 
we must calculate other higher-order diagrams, which are the charm penguin, the vertex 
corrections, and so on. In this study, we consider only the chromo-magnetic penguin, as this 
operator is most sensitive to new physics. 

The decay width for B — > <pK decays is written as 

r = -^^\A\ 2 , (3-i) 

where A is the sum of the leading amplitude and the magnetic penguin amplitude. The 
decay amplitudes from the magnetic penguin are given by 

A MP = V t * b V ts (Mf p + Mt IP + Mf p + Mf p + Mf p + Mf p 

+M? P + Mf? p + Mf 1 p + Mf p ) , (3-2) 

where the indices denote the labels in Fig. The charged modes and neutral modes in 
B — > cf)K have the same contributions of the chromo-magnetic penguin. The amplitudes 
M. are expressed as the convolutions of meson wave functions, a hard part, and the Wilson 
coefficients.^ The meson wave functions are non-perturbative, and the hard part, which 
includes the exchange of a hard gluon, is perturbative. We introduce small transverse mo- 
menta for quarks and anti-quarks in mesons. Large double logarithms are then generated 
by corrections in the meson wave functions. Their resummation leads to the Sudakov fac- 
tor.^ The Sudakov factor guarantees a perturbative calculation of the hard partP^ The 
other double logarithms appear from the end-point region of the parton momenta. The 
resummation of these logarithms leads to the threshold factor in the hard partP^ 

We consider the B meson to be at rest. In the light-cone coordinates, the B meson 
momentum Pi, the K meson momentum P2, and the meson momentum P3 are taken to 
be 

M B , x M B/ , x M B , 9 N , s 

P 1 = -|(1,1,0 T ), P 2 =-^(l-r;,0,Or), P 3 = ^|(rJ,l,0 T ), (3-3) 

where = M^/Mb, and the K meson mass is ignored. The momentum of the spectator 
quark in the B meson is written as k\. The light quark in the B meson has a minus component 
ki given by k± = a^Pf, where X\ is the momentum fraction. The quarks in the K meson 
have the plus components X2P2 and (1 — X2)P 2 + and the small transverse components 
and —k2T, respectively. The quarks in the <fi meson have the minus components x^P^ and 
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(1 — x 3 )P 3 and the small transverse components k 3 T and —k 3 T, respectively. The <j) meson 
longitudinal polarization vector and two transverse polarization vectors e^x are given by 
e = (l/V27>)(-rJ, 1, Or) and e ST = (0, 0, 1 T ). 

For example, the amplitude for Fig. E^a) is written as 

Mf p = -8M B C 2 F ^^ / dxidx 2 dx 3 [ hdb^dhhdhM^iM) 

Jo JO 

x [-(1 - x 2 ) {2cj)i(x 2 ) + r K (3^(x 2 ) + f K (x 2 )) 
+r K x 2 (<Pk( x 2) ~ 4> T k{ x 2))} <P<p(x3) 
-r>(l + x 2 )x 3 <t>£(x 2 ) (3^(z 3 ) - <f>l(x 3 )) 
-r 6 r K (l - x 2 ) (4>k(x 2 ) ~ 4>k(x 2 )) (3<^(x 3 ) + <^(x 3 )) 
-r A r K x 3 (l - 2x 2 ) {(j) P K {x 2 ) + <j> T K {x 2 )) (3^(z 3 ) - $(z 3 ))] 

x£,(t)iV,{x 2 (l - x 2 )} c hf P (A B, C, h, b 2 , b 3 ) , (3-4) 

where r K = m 0K /M B and m K is the chiral factor defined as m 0K = M\l{m& + m s ). The 
meson distribution amplitudes, and so on, are given in the Appendix iBl ^{22(1 — 

x 2 )}° is the threshold factor, where N t = 1.775 and c = 0.3.^ The evolution factors E g (t) 
are defined by E g (t) = {a s (t)} 2 C 8 G(£) exp[— Ssit) — Sx(t) — S in which Ssit), Sxit) 
and S^(t) are the Sudakov factors given in the Appendix El The hard functions are given 
by 

hf P (A, B, C, h, b 2 , b 3 ) = -K (Bh) K (Cb 3 ) 
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x / d0tan0J o (A&itan0) J o (A& 2 tan0) J (Ab 3 tan9) , (3-5) 



o 



where A, B and C are given on Table ITT1 The values A 2 , B 2 and C 2 are the squares of the 
virtual quark and virtual gluons momenta. We assume that the hard scale t is defined as 



t = max(v^n VWl l/K 1/&2, 1/63) • (3-6) 

In the above definition of t, we have chosen the same definition for t as we employed for com- 
putation of the leading diagrams. There is a slight dependence on t in the chromo- magnetic 
penguin amplitude, as we see in the next section. The expressions for other diagrams are 
summarized in the Appendix O 

§4. Numerical result 

The parameters that we used in this calculation followsP M B = 5.28 GeV, 

M K = 0.49 GeV, M s = 1.02 GeV, m b = 4.8 GeV, m t = 174.3 GeV, f B = 190 MeV, f K = 
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160 MeV, U = 237 MeV, f} = 220 MeV, t b o = 1.54 x 10" 12 sec, r B ± = 1.67 x 10~ 12 sec, 
yl^ D = 0.250 GeV. In addition, we used the value of the chiral factor m 0K = 1.70 GeV, 
and the parameter values in wave functions are as given in the Appendix iBl 

The numerical results are listed in Table HI The leading amplitude for B° —>■ <f)K° is 

Table I. Magnetic penguin contributions for B — > cfiK. 



A4f p + Mf p -0.0308 - i 0.0169 

Mf p + Mf p 0.00309 + i 0.00544 
Mf p + Mf p negligible 
+ Mf p + M™ p 

Mf p -0.00229 -i 0.000798 

Mf p 0.00623 -i 0.0147 

M MP -0.0238 -i 0.0270 



0.0944 -i 0.0383. 181 We find that the diagrams (a) and (b) are dominant, and the magnetic 
penguin amplitudes are about 30% of the leading amplitudes. 

We study the hard scale dependence of chromo-magnetic amplitudes. We define the hard 
scale as t = k ■ max(y / pP], ^J\B 2 \, 1/bi, l/b 2 , 1/&3), where k is a parameter, and vary k 
between 1 and 1.5. Figure 121 shows the hard scale dependence of -M* /p + M.^ IP . We find 
that there is a slight dependence on t in the chromo-magnetic penguin. This is consistent 
with the fact that the Wilson coefficient for 0$g is most sensitive to new physics, and its 
contribution, by itself, is physical. 
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Fig. 2. t dependence of M^ IP +M^ IP . The 
lower curve is the real part and the upper 
curve is the imaginary part. We see that 
there is a slight dependence on t in the 
chromo-magnetic penguin. 




0.4 0.6 

Fig. 3. The distribution of q 2 for ReM™ p . 
The expectation value of q 2 is 6.3 GeV 2 . 



The value of q 2 of the virtual gluon in 0' 8G is written in terms of x\, £2 and X3, which are 
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the momentum fractions of partons. For example, in the diagram (a), q 2 = (1 — X2)x^M\. 
In Fig. El we show the distribution of q 2 for Ai^f . The shape of this graph is not simple. 
The expectation value of q 2 is calculated as 

(q 2 ) = 6.3 GeV 2 . (4-1) 

This is approximately equal to M§/4. 

Next, we analyze the FA using the value of q 2 given in Eq. (j4-lj) . In the FA, the effective 
Hamiltonian for the chromo-magnetic penguin is^ 

We use the value 6.3 GeV 2 for (q 2 ) and 1.5 GeV for the renormalization scale, which is 
a typical scale of the hard part in the PQCD approach. The result for the amplitude 
is M. = —0.029, where we input the value Fbk — 0.38. This is to be compared with 
m mp + m mp = _q.0308 - i 0.0169. The real part is agree with the above result. 

§5. Summary 

In this paper, we calculated the chromo-magnetic penguin amplitudes for B — > cf)K decays 
in the PQCD approach. We find that the chromo-magnetic penguin contribution is about 
30% of the leading contribution. Therefore, we find that the chromo-magnetic penguin is 
not negligible in hadronic two-body decays. However, for the Standard Model there are 
many other higher-order diagrams that must be considered simultaneously if we are to add 
the magnetic penguin term. If we consider the contribution of new physics, this diagram 
can be considered separately. It is noted that there is only weak energy scale dependence. 
In addition, we calculated the expectation value of q 2 of the gluon in the chromo-magnetic 
penguin, and the result is approximately M B /4. 
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Appendix A 

The Sudakov Factor 



The Sudakov factor is written aa*' 



S B {t) = s{x 1 P^,b 1 ) + / — l^sip)) , 

Ji/h V 

S K (t) = s(x 2 P 2 +, 62) + s((l - a? 2 )P+, 62) + / —^(a s (fi)) , 

= s(x 3 P 3 -, 63) + s((l - x 3 )P 3 ~, 63) + / —^(a s (fx)) . 

Jl/b 3 l 1 



The exponent s is given by 



ln( ^ ) A(a.0x)) + B(a.0x)) 



where the anomalous dimensions A, to two loops, and B, to one loop, are 



7T 



9 3 27 3 V 2 



7T 



5 = -— In ' 

3 7T V 2 



2a s , ( 



(Al) 
(A-2) 
(A3) 



(A-4) 



(A5) 
(A6) 



with Cf = 4/3 being a color factor, / being the number of active flavors, je being the Euler 
constant, and (3q = (33 — 2/)/3. The anomalous dimension of mesons is given by 



7*(«-(/*)) = 2 %M^)) = ~ 

Appendix B 

Wave Functions 



7T 



(A7) 



The B meson wave function is defined by 



dxd 2 k T e- l(xP ~ z+ - kTZ ^ [( f> + M B ) l5 cj> B {x : k T )} 



aj3 ' 



(B-l) 



where the indices a and (3 are spin indices, % and j are color indices, and N c is the color 
factor. The distribution amplitude 0b is normalized as 



B 



dxrfuixiM = 0) = 

Zy/Zl\ c 



(B>2) 
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where 61 is the conjugate space of hi, and fs is the decay constant of the B meson. In this 
study, we used the model functions 



!>s(x, b) = Nbx 2 (1 — x) 2 exp 



xM 



B 



u\b 2 



(B-3) 



where Nb is the normalization constant and uob is the shape parameter. We used ujb = 0.4 
GeV, which was determined by the calculation of form factors.^ 
The K meson wave function is given by 



-iSij r 1 



dxe ixP ^ 5 [P<I>k(x) +m m (j) P K (x) +m 0K (/j jh - l)(j) T K {x)\ 



;b-4) 



where m K = M^/i^ma + m s ), n M = z^/z~~ and = P^/P + . The function <j*j({x) is the 
leading-twist contribution and 4>k an d $k are twist-3. These distribution amplitudes were 
calculated using the light-cone QCD sum rule™ 



<f> A K {x) 



f. 



X 



2^/2N c 
fx 



4> T K {x) 



2^2N C 
fx 



2V2N, 

x 



6x{l - x) 1 + aiCf(l - 2x) + a 2 Cf (1 - 2x) 
l+Uor ]3 -^p 2 K \cl(l-2x) 

-3 |r/3^3 + + 6a 2 )| C/(l - 2x) 

1 - 2x) 



(B-5) 



;b-6) 



1 ( 5^ - ^3^3 - ^Px - jjp^a 2 ) ( 1 - 10./' 1 0.r- 



(B-7) 



Here, = (m^ + m s )/M K and C^(x) is the Gegenbauer polynomial. The parameter values 
are ai = 0.17, a 2 = 0.20, 773 = 0.015 and u 3 = —3.0. 

The cf) meson wave function of the longitudinal parts is given by 



5. 



dx^e 



[M+ 4^M x s)+ U M+M + m^%{x z )] 



aj3 



(B-8) 



Here, <^ is the leading-twist distribution amplitude, and the other terms are the twist-3 
amplitudes. They were also calculated using the light-cone QCD sum ruleP^ 



u 

2y/2N ( 



:6x(l — X) 



;b- 9 ) 
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f q 



2^/2N c 



■Is 



3(1 - 2a:) 2 + yC 3 T {3 - 30(1 - 2xf + 35(1 - 2xf} 
+^ + |l-(l-2a;)log^ 

(1 - 2a?) {6 + 95+ + 140C 3 T (1 - 10x + 10x 2 )} 

x 



-35+ log 



I — x 



Here, (I = 0.024 and 5+ = 0.46. 



(B-10) 



(B-ll) 



Appendix C 

Magnetic Penguin Amplitudes — 

Below, we present the decay amplitudes for Fig. |H(b) - (j): 



Mff p = 8M 6 B C 2 F - 



M 



MP 



MP 



2N, 



c JO 



dx±dx 2 dx 3 / bidbib 2 db 2 b 3 db 3 (f)B(xi,bi 



(CM) 



x [x 1 (t)K(x 2 )(t)4 > {xz) - 2r K (-2 + Xi)^^)^^) 
+r 4> x 1 <t>i(x 2 ) (30J(ar 3 ) + $(a*)) 

-12r r K xi0^(x 2 )^(a;3) + 2r K x 3 <f K {x 2 ) (30£(ac 3 ) - 4>\{x 3 ))] 

xE g (t)N t { Xl (l - x 1 )} c hf p (A, B, C, b 2 , h, b 3 ) , 

(N \ y/2N f°° 
C F - -y j y N ° J dx 1 dx 2 dx 3 J b l db 1 b 2 db 2 b 3 db 3 4> B (x 1 ,b 1 ) 

X [2(1 - X 2 )(l -Xx- X 3 )0^(X2)0^(X 3 ) 

+r K (l - x 2 )(l -X1-X2- x 3 ) {(p P K {x 2 ) - 4> T K {x 2 )) (pcj,(x 3 ) 
+7>a? 3 (l - ari - x 2 - x 3 )(f)^(x 2 ) (0J(ar 3 ) + ^(x 3 )) 
+r#r> {(1 - xi)(l - x 2 ) - x 3 } (0^(x 2 ) - <p T K {x 2 )) (3<^(x 3 ) + <^(x 3 )) 
+2r K r (j) x 2 x 3 (3^(x 2 )^(x 3 ) - 0^(a; 2 )^(x 3 ))] 
xE 5 (t)C P (A, B,C, b u b 2 ,h), (C-2) 

b 1 db 1 b 2 db 2 b 3 db 3 (f) B (x 1 , b x ) 



c ao 



dx\dx 2 dx z 



x(xi -x 3 ) [-(2 + x 2 )0^(x 2 )0 (x 3 ) 

{ (0/^(^2) - 0X^2)) - X 2 (0^(X 2 ) + 0x(x 2 )) } <MX 3 ) 

+r> (2xi - x 2 - 2a; 3 ) 0^(x 2 ) (<^(x 3 ) - <^(x 3 )) 

+2r*7> { (0£(x 2 ) - 0^(x 2 )) - x 2 {(j) P K (x 2 ) + <$(x 2 )) 4(x 3 )}] 



(C-3) 
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Mf p = 8M|C|^P f dx x dx 2 dx 3 f°° hdbibtdhhdbsfefaM) 

2^c Jo Jo 

x [3x30^(x 2 )0 (x 3 ) + 2r K (2 - x ± + x 3 ) P ;(x 2 )0 (x 3 ) 
+r 4> x 3 {1-2 ( Xl - x 3 )} <f>i(x 2 ) (<t>%{x 3 ) - <j>l(x 3 )) 
-Qr K r 4 ,x 1 (j)K(x2) (tyM - 4>\M) + l^rxr^^M^M] 

xE g (t)N t {x 3 (l - x 3 )} c h™ P (A, B, C, 62, 63, 61) , (C-4) 

Mf p = -&M%C 2 F ^E± [ d Xl dx 2 dx 3 f°° b 1 db 1 b 2 db 2 b 3 db 3( j> B (x 1 ,b 1 ) 
2N C J J 

x [-3xi(l - x 2 )<j>K(x2)<j>4,(x 3 ) 

+r K (l - x 2 )(2x 1 - x 2 ) ((Pk(x 2 ) + 4> T K {x 2 )) 4>4,{x 3 ) 

+2r <t> {x l - 2x 2 )0£(x 2 )0^(x 3 ) 

-Qr K r^ Xl (l - x 2 ) (0£(x 2 ) - (t> T K {x 2 )) 0^(x 3 ) 

+Qr K r <p x 2 (0£(x 2 ) + 4> T K {x 2 )) <j>%{x 3 )\ 
xE g (t)N t {x 2 (l - x 2 )} c h™ p (A, B, C, b 3 , b 2 , h) , (C-5) 

Mf p = 8M B C F (c F - jf d Xl dx 2 dx 3 jf hdhhdhhdhMxi, 61) 

x (xi - x 3 ) [-(1 - x 2 )0^(x 2 )00(x 3 ) - (^(^2) - 0k (^2)) 00(2:3) 
-2r^0^(x 2 )0 (a;3) + r>(l + x 2 )0^(x 2 ) (0^(x 3 ) - 0^ 3 )) 
-4r 0£(x 2 )0^(x 3 ) 

+2r K r (j)X2 (0£(x 2 ) - 4>k(x 2 )) 0* (x 3 ) + 4r x r </ ,0 T -(x 2 )0^(x 3 )] 
xE a (t)h*f p (A, B,C,b 2 , 63, 61), (C-6) 
p = 8M|C F f C F - ^ J d Xl dx 2 dx 3 J™ hdb^db^dhMxi, 61) 

x [-zi(l - a; 2 ) 0^(^)00(2:3) - ^^(l - x 2 ) (0^^) + 4> T K {x 2 )) 0</,(x 3 ) 
-r x x 2 (2 - xi) (<Pk(x 2 ) - <Pk(x 2 )) (x 3 ) 
+r^xi(xi + x 3 )0 x (x 2 ) (0^(x 3 ) - 0* (x 3 )) 

+7> {-X! + 2X 2 (1 - X! - X 3 )} 0^(X 2 ) (0^(X3) + 0;(X 3 )) 

-r x r^x 2 (l - x3) (0£(x 2 ) + 0^(x 2 )) (30* (x 3 ) + 0' (x 3 )) 
-2r^r x 2 (0 P : (x 2 ) + x (x 2 )) 0^(x 3 ) 
+2r x r xix 2 (30 P : (x 2 )00(x 3 ) + 0x(x 2 )0^,(x 3 )) 
-2r K r^ Xl (0 P : (x 2 ) - 0^(x 2 )) <^(x 3 )] 
*E a {t)h% p (A,B,CMMM) , (C-7) 
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M? IP = -4M 4 B ^^C F N C [ dx lC lx 2 dx 3 [ b 1 db 1 b2db 2 b 3 db 3 <f>B(xiM) 

-r$4>i(x 2 ) (<f>%{x 3 ) + 4>l(x 3 )) +Ar^ K <p p K (x 2 )<Pl{xz)\ 
xE g (t)/ i f p ( J B,C,6 1 ,6 2 ,6 3 ) , (C-8) 

Mf p = 2M%^^C F N C [ d Xl dx 2 dx 3 [ b 1 db 1 b 2 db 2 b z dh(j) B {x 1 M) 
27V C J J 

x [-2 {xix 2 - (2 + x 2 )x 3 } 0i(x 2 )00(x 3 ) 

+r K {3(2 + xi(-2 + x 2 ) + x 3 - x 2 x 3 )0^(x 2 ) 

+ (2 + xi(2 + x 2 ) - Zx 3 - x 2 (4 + x 3 ))4> T K (x 2 )} ^(x 3 ) 

+r (f) {xi - 2xi(-3 + x 2 + 3x 3 ) + x 3 (-2 + 2x 2 + 5x 3 )} </4(x 2 )<^(x 3 ) 

-7>(-2 + 3xi - 2x 2 - 3x 3 )(x 1 - x 3 )(t)j ( {x 2 )(j) t 4) {x 3 ) 

+6r r x {3(-Xi + x 3 )0^(x 2 ) - (xi - xix 2 + x 2 x 3 )0^(x 2 )} <^(x 3 ) 

-2r^r K {-3(a;i - Xia; 2 + x 2 x 3 )0^(x 2 ) 

+ (xx(3 + 2x 2 ) - (1 + 2x 2 )x 3 )0^(x 2 )} <f>l(x 3 j\ 
xE g (t)hf p {A,B,CMMM) • (C-9) 

The hard scale t is defined in Eq. ()3-6j) . and the hard function h^ IP is defined in Eq. ()3-5|) . 
The other hard functions are defined as follows: 

^f P (A B, C, &X.&2, fog) 



b x db x b 2 db 2 b 3 db zi —^ / dd\dd' % dd^ 



xifo ^C^/6? + 6i + 26 1 6 2 cos(0' 1 -^ K (b^I$ + &§ + 26 2 6 3 cos(£ 2 - 9> 3 ^j 



^XcosW^+^sinP^ (C10) 







hf p (B } CMMA 



3/ 



1/(2^) for &1 + 6 2 >6 3 >|6 1 _6 2 | 
K {Bh) K {Cb 3 ) < . , (C-ll) 

I otherwise 

hf p {A,B,C } b x ,b 2 M) 

,2. / & 2 , v /6? + ^ + 26 1 6 2 cos(^-0 2 )) 

^ 

xKo^&OiTo f + &| + 2b 1 b 2 cos(9 1 - 6 2 )J . (C-12) 
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Here, A = \Zlb\h\ + 26p| + 2b\b\ - 6f - b\ - &J/4 and D = b\ cos6' x + 6 2 cos0£ + 63 cos# 3 . 
The quantities A, £? and C, which are functions of x\, X2 and x 3 , are given in Table ITT1 For 
simplicity, we ignored \kix — k 3 x\ in the calculation of hf Ip . 



Table II. The definitions of A, B and C 



Diagram 


A 


B 


C 


(a) 


y/x~ 2 M B 


^XiX 2 M B 


Za/(1 - x 2 )x 3 M B 


(b) 




^/xiX 2 M B 


a/xi - x 3 M B 


(c) 


i\/x 2 (l -x\- x 3 )M B 


yJx x x 2 M B 


za/(1 - x 2 )x 3 M B 


(d) 


\fx 2 {xi - x 3 )M B 


sJx\X 2 M B 


\]x\ - x 3 M B 


(e) 


i^fx~ 3 M B 


-x 2 )x 3 M B 


sjx x - x 3 M B 


(f) 


- x 2 M B 


-x 2 )x 3 M B 


^x x x 2 M B 


(g) 


y/(xi - X 3 )(l - X 2 )M B 


?a/(1 -x 2 )x 3 M b 


sjx\ - x 3 M B 


(h) 


a/1 - x 2 (l -x x - x 3 )M B 


i<s/(l -x 2 )x 3 M B 


y/XlX 2 M B 


(i) 




^Jx x x 2 M B 


Za/(1 -x 2 )x 3 M b 


0) 


y/XlX 2 M B 


2a/(1 -x 2 )x 3 M b 


y/Xi - X 3 M B 
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